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ABSTRACT 



In a hierarchical cosmogony galaxies build up by continuous merging of smaller 

structures. At redshift three the matter content of a typical present-day galaxy is 
dispersed over several individual clumps embedded in sheet-like structures, often 
aligned along filaments. We have used hydrodynamical simulations to investigate the 
spatial distribution and absorption properties of metal enriched gas in such regions of 
ongoing galaxy formation. The metal and hydrogen absorption features produced by 
the collapsing structures closely resemble observed QSO absorption systems over a wide 
range in Hi column density. Strong Cll and SilV absorption occurs for lines-of-sight 
passing the densest regions close to the center of the protogalactic clumps, while CiV is 
a good tracer of the prominent filamentary structures and Ovi becomes the strongest 
absorption feature for lines-of-sight passing through low-density regions far away from 
fully-collapsed objects. The observed column density ratios of the different ionic species 
can be well reproduced if a mean metallicity [Z/H] =—2.5, relative abundances as found 
in metal-poor stars, a UV background with intensity J-22 = 3 at the Lyman limit, and 
either a power law spectrum (J oc v^^'^) or the spectral shape proposed by Haardt 
&: Madau (1996) are assumed. The observed scatter in [C/H] is about a magnitude 
larger than that in the simulations suggesting an inhomogeneous metal distribution. 
Observed and simulated Doppler parameter distributions of Hi and Civ absorption lines 
are in good agreement indicating that shock heating due to gravitational collapse is a 
second important heating agent in addition to photoionization heating. The two-point 
correlation function of Civ lines agrees reasonably well with the observed correlation 
function suggesting that its high-velocity tail is caused by occasional alignments of 
several absorbing clumps still expanding with the Hubble flow. Both high-ionization 
multi-component heavy-element absorbers and damped Lyman alpha systems can 
arise from groups of moderately sized protogalactic clumps (Mbaryon~ W^Mq). Recent 
detections of star forming galaxies at similar redshifts are consistent with this picture. 
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1. Introduction 



Can we understand galaxy formation at high redshift studying QSO absorption systems? 
Lines-of-sight (LOS) to distant QSOs must frequently traverse protogalactic environments. If 
the gas in these regions has already experienced some metal enrichment the galaxy formation 
process should produce a characteristic pattern of metal absorption lines in the QSO spectrum. If 
we succeed in identifying the spectroscopic signature of such regions, high-resolution spectra of 
QSOs at high redshift should be able to yield detailed insights into the physical state of the gas in 
young galaxies and the gaseous reservoir from which they form. Once the observable absorption 
features have been related to the properties of the absorbing objects one may hope to be able 
to actually constrain the galaxy formation process and cosmological parameters. This paper is 
mostly concerned with the first aspect: what are absorption systems (in particular: heavy element 
systems) at high redshift, and can we model them realistically by adopting a specific scenario of 
structure formation? 

QSO absorption studies provide only one-dimensional information, so observations of single 
LOS cannot give independently the spatial density and the size of the objects causing the 
absorption. Only the product of number density and cross-section is constrained. Additional 

information about the spatial distribution of the absorbing gas has to be sought, e.g. by the 
identification of heavy element absorbers with a class of galaxies whose number density is known. 
As Bahcall & Spitzer (1969) and Burbidge et al. (1977) pointed out, the cross section of the 
absorbing gas would have to be much larger than typical half-light radii of present-day spiral 
galaxies if these objects had the same comoving space density as such galaxies. First observational 
evidence in favor of large metal absorption cross-sections came from work by Bergeron (1986) who 
discovered low redshift galaxies close to the LOS to QSOs with known Mgii absorption systems, 
coincident in redshift with the absorbers. Surveys of optically thick (Lyman Limit) absorption 
systems and galaxies at intermediate redshift have reported results consistent with galaxies being 
surrounded by Mgii absorbing halos of radius ~ 40 kpc (e.g. Bergeron (1995), Steidel (1995), 
Churchill, Steidel Sz Vogt (1996)). Extending such work to optically thin systems and lower 
redshifts Lanzetta et al. (1995) found evidence for even larger Hi halos with radii of order 160 
kpc. Similarly, damped Lya systems were interpreted by Wolfe and collaborators (e.g. Wolfe 
1988, Lanzetta ct al. 1991) as large protodisks, progenitors of present-day spiral galaxies with 
significantly larger radii at high redshift. 

The large-halo/disk scenario can qualitatively explain the component and internal structure 
of heavy element absorption systems, especially the strong clustering measured for CiV systems 
(Sargent et al. 1979, Sargent, Boksenberg & Steidel 1988, Petitjean Sz Bergeron 1994). In this 
picture the individual absorption components could be clouds orbiting in a halo or co-rotating in 
a disk, produced and replenished e.g. by thermal instability (Bahcall 1975, Fall & Rees 1985, Mo 
1994, Mo & Miralda-Escude 1996). 

However, unambiguous evidence in favor of large absorption cross-sections and the 
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identification of the absorbing objects with massive galaxies is restricted to low redshift 
observations. Little is known about the sizes of damped Lyct absorbers at high redshift (e.g. 
M0llcr & Warren 1995, Djorgovski et al. 1996) and massive disks arc not the only viable dynamical 
model for the observed velocity structure. Moreover, observations of galaxies at high redshift 
are strongly biased towards objects with high current star formation rates. Thus it is not clear 
whether the transverse separations on the sky of galaxy-absorber pairs coincident in redshift do 
reliably indicate the presence and sizes of any (hypothetical) halos/disks. It is difficult to ascertain 
that there is not an undetected (not necessarily faint) galaxy closer to the LOS. 

The possibility of "smaller" but more numerous objects as sources of the metal absorption 
has been explored earlier. Tyson (1988) suggested identifying damped Lya systems with gas-rich 
dwarf galaxies instead of large proto-disks. York et al. (1986) discussed clusters of such objects to 
explain the component structure of Civ systems. Individual galaxy halos cannot produce potential 
wells deep enough to explain the largest velocity splittings of Civ systems (up to 1000 km s~^) 
as virialized motions. Pettini et al. (1983) and Morris et al. (1986) concluded that the velocity 
splitting and large cross-sections are equally difficult to understand if objects similar to present-day 
galaxy clusters placed at high redshift were causing the absorption. 

Earlier attempts at understanding QSO absorption systems have mostly employed heuristic 
models. The ionization state of the gas was calculated with simplified assumptions about the 
geometry, temperature and dynamical structure of the gas. Recently, realistic hydro dynamical 
simulations of the fate of gas in a universe subject to hierarchical structure formation (Cen et 
al. 1994, Hernquist et al. 1996) have led to a deeper understanding of the large scale distribution 
of baryons. In the new picture a coherent filamentary large scale structure of dark matter and 
baryons is responsible for Lya forest absorption lines (e.g. Petitjean, Miicket & Kates 1995; 
Zhang, Anninos &; Norman 1995; Miralda-Escude et al. 1996). Denser condensations embedded in 
the filaments are unstable against rapid collapse and cooling of the gas and probably form stars at 
their centers. 

In such a hierarchical picture (cf. Lake 1988) metal absorption systems at high redshift are 
more likely to arise from groups of relatively small, merging protogalactic objects, rather than 
from ensembles of clouds in huge virialized halos: 

(1) Recent Keck spectroscopy of Lya forest systems (Cowie et al. 1995, Tytler et al. 1995) 
has shown that contamination of Lya forest clouds by carbon is common in high-redshift Lya 
absorbers with Hi column densities as low as a A'^(IIi) ~ a few xlO^^cm"^. In the numerical 
simulations such column densities typically correspond to gas densities smaller than those expected 
for fully collapsed objects at these redshifts (indicating baryonic overdensities of order ten), so at 
least some metal absorption systems appear to occur outside virialized regions. 

(2) Typical observed temperatures of Civ systems are somewhat larger than expected if 
the gas were heated only by photoionization (Ranch et al. 1996, hereafter RSWB). Such an 
enhancement is predicted by numerical simulations (Haehnelt, Steinmetz & Ranch 1996, hereafter 
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paper I) and is most likely due to shock and (to a smaller extent) adiabatic heating during the 

gravitational compression of the gas. Absorber sizes along the LOS inferred from ionization 
balance calculations are of order ten kpc. This is uncomfortably large for a cloudlet-in-halo model 
(Hachnclt, Ranch, & Stcinmetz 1996, hereafter paper II). 

(3) Simple scaling laws predict that in a hierarchical scenario the typical ratio of cooling time 
to dynamical time decreases as (1 + z)~^/^. At high rcdshift the gas will generally cool rapidly 
out to the virial radius. It will be difficult to maintain large, massive, hot halos for an extended 
period of time. Similarly for fixed circular velocity both the mass and the size of the dark matter 
component of typical objects at high redshift decrease as (1 + z)"^/^ (see Kauffmann (1996) and 
Steinmetz (1996a) for quantitative predictions). 

In paper I we used numerical hydrodynamical simulations to demonstrate that observed 
Civ absorption systems can be qualitatively understood in terms of absorption by protogalactic 
clumps (PGCs) formed through gravitational instability in a hierarchical cosmogony. The Civ 
and Hi absorption features caused by groups of PGCs were found to resemble observed Civ and 
Hi absorption systems if an overall homogeneous abundance of [C/H] —3 to —2 was assumed. 
Here, we will perform a quantitative analysis and extend the work to a larger set of diagnostically 
useful ionic species. Section 2 will give the details of the numerical modelling. In section 3 we 
discuss the line formation process and the physical nature of the absorbing structures in regions of 
ongoing galaxy formation. In section 4 we investigate the observational consequences. Conclusions 
are drawn in section 5. 



2. Numerical simulations 
2.1. The code 

The simulations were performed using GRAPESPH (Steinmetz 1996b). This code combines 
a direct summation N-body integrator with the Smoothed Particle Hydrodynamics technique 
(Lucy 1977, Gingold &; Monaghan 1977) and is implemented for the special purpose hardware 
GRAPE (Sugimoto et al. 1990). The version of the code we use here is especially adapted to 
follow a mixture of collisionless (dark matter) and collisional (gas) fluids. It is fully Lagrangian, 
three-dimensional, and highly adaptive in space and time as the result of the use of individual 
smoothing lengths and individual timesteps. The physical processes included in this version of 
the code include self-gravity, pressure gradients, hydrodynamical shocks, radiative and Compton 
cooling and photoheating by a UV background with a specified spectrum. We do not assume 
rate equilibrium for the dominant baryonic species (H, H+, He, He+, He^^, and e~) but follow 
self-consistently their non-equilibrium time evolution. However, this had only modest effect on the 
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gas temperature in the lowest density regions. The non-equihbrium time evolution can generally 
be neglected for the questions addressed in this paper but will be important for low column density 
systems at epochs close to rcionization. Wc assume that the gas remains optically thin throughout 
the calculation and that the background radiation is uniform in space. Radiative transfer effects 
have been neglected (see also Navarro & Steinmetz 1996). 



2.2. Initial conditions 

The initial conditions are identical to those in Navarro and Steinmetz (1996) and were 
originally designed to study the formation of galaxies with circular velocities between 80 km s~^ and 
200kms~^. The background cosmogony is a CI = 1, Hq = 50kms~^Mpc~^ cold dark matter 
(CDM) universe with a normalization of as = 0.63. The baryon fraction is Cl^ = 0.05. Based on a 
P3M large scale structure simulation in a 30 Mpc box, eight halos were selected (four each with 
circular velocity of about 100 km s^^ and 200 km s^^) and resimulated at much higher resolution. 
The tidal fields of the large scale matter distribution of the original simulation were still included. 
The high-resolution sphere has a radius of about 2.3 (3.8) Mpc (comoving) for the systems with 
a circular velocity of 100 km/sec (200 km/sec). The mass per gas particle is 4.9 x IO^Mq and 
2.3 X 10^ M0 in the low- and high-mass systems, respectively. We adopt a Plummer gravitational 
softening of 2.5 (5) kpc for the dark matter and of 1.25 (2.5) kpc for the gas particles in the low 
(high) mass systems. All runs are started at 2; = 21. 

We use several different description for the intensity, redshift dependence and frequency 
dependence of the UV background. Most of the simulations were done using a power-law spectral 
energy distribution and the redshift dependence suggested by Vedel, Hellsten & Sommer-Larsen 
(1994). We varied turn-on redshift, spectral index (a = 1 . . . 2) and normalization of the 
background {J-22 = 0.3 . . . 30, in units of 10~^^ ergs~^ cm"^ Hz~^ sr~^). The more realistic 
description of the UV background proposed by Haardt & Madau (1996) was also used. Although 
varying the UV background had noticeable effects on the temperature and the ionization of the 
gas, the effects on the overall gas distribution were small. In this paper we will only briefly 
discuss these effects and concentrate on simulations with a = 1.5 and J-22 = 3 for which the 
background was switched on at the beginning of the simulation {z = 21). This corresponds to a 
typical quasar-like spectrum and is consistent with measurements of the helium Gunn-Peterson 
decrement (Davidsen, Kriss &; Zeng 1996). The initial gas temperature was assumed to take the 
value, for which photoionization heating balances the cooling due to collisional processes and 
Compton scattering. For simplicity, wc take the gas to be homogeneously contaminated by metals 
(C,Si,N,0) with solar relative abundances and an absolute metallicity of 10~^ solar. Effects of time 
dependent and/or inhomogeneous chemical enrichment were neglected. Later we shall consider 
different total and relative metallicities and briefly discuss implications of an inhomogeneous metal 
contamination. 
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Fig. 1. — The upper panel shows a grey-scale plot of the projected Hi column density (logA^) of 
the inner 700 kpc of a simulation box at z = 3.07 which will contain three Vc ~ 100 kms~^galaxies 
at 2; = 0. The lower panel shows the temperature (log T, Hi column density weighted) for the same 
box. 

2.3. The matter and temperature distribution 

Some of the simulations were run up to the present epoch. For an individual set of initial 
conditions, typically between one and three disk like galaxies with circular velocities between 80 
and 200kms~^have formed by now. Each consists of several thousand particles. In addition, a 
couple of halos with smaller masses and circular velocities (about 50kms^^) are present. Size, 
mass and circular velocities of these disks are in good agreement with observed spiral galaxies. 
Compared to present-day spiral galaxies, however, the disks are too concentrated. This may be an 
artifact due to the neglect of energy and momentum input by star formation and active galactic 
nuclei. 

Tracing the evolution of a present-day galaxy backwards in time, one finds that at redshift 
three it has "split" into several progenitors separated by several hundred kpc. Generally the 
circular velocities of these progenitors are slightly smaller than that of the galaxy which they will 
later form, but this is not true for the most massive progenitors. These have the same or even 
higher circular velocities even at redshifts of three to five, although their mass may be a factor ten 
smaller (Steinmetz 1996a). The progenitors probably suffer from the same overcooling problem 
as their present-day counterparts. This will have small effects for intermediate column density 
absorbers but conclusions regarding damped Lya absorbers should be drawn with some caution. 
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Figure || shows the column density and temperature distribution (Hi weighted along the LOS) 
at z = 3.07 (bright=high column density/hot) of a typical simulation box in which three galaxies 
with circular velocity Vc ~ 100 km will have formed at redshift zero. The area shown is 700 kpc 
across (proper length). Individual PGCs are embedded in a filamentary network. Regions leading 
to a single Vc ~ 200kms~^ galaxy look similar, with individual components scattered over several 
hundred kpc. 

The temperature in the gas rises from less than 10^ K at the center of expansion-cooled voids 
(large diffuse dark areas) to a few times 10^ K for the gas in the filaments. Higher temperatures 
up to several 10^ K degrees occur in spheroidal envelope within ~ 30 kpc of the PGCs. These 
envelopes of hot gas arise when infalling gas is shock-heated to temperatures of the order of the 
virial temperature, but the density of this gas is low enough that cooling times are long compared 
to dynamical time scales. The temperature, the location and the mass fraction of this hot gas 
component can depend sensitively on the assumed spectral shape of the UV background and on 
the metal enrichment (Navarro & Steinmetz 1996). In the innermost few kpc of the PGCs the 
cooling timescales are always short compared to the dynamical time scale and the dense gas cools 
precipitously to temperatures below 10'* K. 



2.4. The ionization state of the gas 



We have used the photoionization code CLOUDY (Ferland 1993) to calculate the ionization 
state of the gas for the temperature, density, and UV background of each SPH particle in the 
simulation. We assumed an infinite slab of gas of low metallicity, optically thin to ionizing 
radiation and illuminated from both sides by a homogeneous UV field. (Self-)shielding from UV 
radiation in optically thick regions was generally not taken into account. Thus quantitative results 
for LOS with column density above 10^ cm~^ need to be viewed with some caution. Below we will 
concentrate on the ionic species Hi, Cii, Civ, Siiv, Nv, Ovi which are known to show the strongest 
observable lines in low and intermediate column density absorbers. As discussed in papers I and 
II, observations and SPH simulations indicate that gas temperatures can significantly deviate from 
the equilibrium value where photo-heating balances line cooling processes. It is then important to 
use the actual gas temperature to calculate the ionization state of the different states. Figure |2| 
shows the ratios of various metal to Hi column densities as a function of density for our chosen 
ionic species and a set of temperatures between 10* and 10^ K . At high and intermediate densities 
the metal column density generally increases relative to the Hi column density with increasing 
temperature while at small densities the opposite is true. Siiv and Civ are especially temperature 
sensitive in the temperature and density range most relevant for intermediate column density 
absorption systems around njj ~ 10~* cm~^ and T ~ 5 x 10* K. 
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Fig. 2. — Column density ratios relative to Hi for six different ionic species calculated with 
CLOUDY for fixed density and fixed temperatures between 10"^ K and 10® K as indicated on the 
plot. A power law with a = —1.5 and J-22 = 3 was assumed for the background UV field. 
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Fig. 3. — Projected column density of the inner 700 kpc of a simulation box at z = 3.07 which will 
contain three Vc ~ 100kms~^galaxies at z = 0. Shown are logarithmic column density contours in 
steps of 1 dex for Hi, Cii, Civ, Siiv, Nv and Ovi. [Z/H] = -2 and solar relative abundances were 
assumed. 
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Figure ^ shows logarithmic column density contours for all 6 ions in steps of 1 dex. Cii 
and Siiv are confined to dense regions marking likely sites of star formation in the center of the 
PGCs. Civ traces the hotter filaments and halos, while Ovi matches well the low column density 
Hi contours and traces low density gas of any temperature (cf. Chaffee et al. 1985, 1986). In the 
box shown the covering factor for detectable hydrogen absorption is 100%. The same is true for 
Ovi if the metal contamination is indeed homogeneous and does not drop in regions outside the 
filamentary structure. Approximately one third of all LOS give detectable Civ absorption. 



2.5. Generating artificial QSO spectra 

To analyze the appearance of the galaxy forming region in absorption we drew 1000 LOS and 
generated artificial spectra for the simulation box shown in figure 1 (see section 2.3). The LOS 
had random orientations and random offsets within ±225 kpc from the center of the box. Optical 
depth t{v) profiles along the LOS were constructed projecting the Voigt absorption line profile 
caused by the column density of each individual spatial pixel onto its proper position in velocity 
space, using the relation (for the Hi Lya line) 

riv) = y^n{v - V,) = 1.34 X 10-12 i^(^;) (1) 

(cf. Gunn &: Peterson 1965). Here Tiiy — Vi) is the optical depth of the Voigt profile at the 
observed velocity v, caused by the column density AA'^j in spatial pixel i, moving with velocity 
Vi = Vpec{i) + vuuhhieii) ■ This can be expressed in terms of the column density per unit velocity, 
dN/ dv at V. Units of N are cm"^, v is in kms"^. 

The spectra were made to resemble typical Keck data obtainable within a few hours from a 
16-17th magnitude QSO (S/N=50 per 0.043 A pixel, FWHM = 8km/s). These "data" were then 
treated and analyzed in exactly the same way (by fitting Voigt profiles) as actual observations. An 
absorption line was deemed detectable when the equivalent width had a probability of less than 
IQ-^ (corresponding to a > 4.75 a event) to have been caused by a statistical fluctuation. 



3. Absorption line formation in regions of ongoing gravitational collapse 

Having access to the three-dimensional gas/temperature distribution and the velocity field of 
the gas we are able to study the line formation process as a function of the characteristic physical 
properties of the absorbing structures. We expect the general mechanisms to be similar to those 
found earlier for Hi absorption features (Cen et al. 1994; Zhang et al. 1995; Hernquist et al. 1996; 
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Miralda-Escude et al. 1996). The emphasis of the present work is, however, on the much narrower 

absorption lines from metal ions. These allow us to study the velocity structure in galactic 
potential wells, where the corresponding Hi Lya lines are saturated and strongly blended. The 
resolution of the present simulations is much higher than that of previous work, so we can pursue 
the physical quantities into regions of larger densities, and study the fate of objects of small mass. 



3.1. Environments causing specific absorption patterns 

The following figures illustrate a few "typical" situations where gravitational collapse gives 
rise to Civ and other metal absorption features. The panels give (from top to bottom) the spectra 
of the six ions offset vertically by 0.5 , the total baryon density, the temperature and the peculiar 
velocity. The x axis is the spatial coordinate along the LOS labeled by the relative velocity each 
position would have if following an undisturbed Hubble flow. In the case of Civ, Siiv, Nv, and 
Ovi only the stronger transition of each doublet is shown. The dotted lines connect the spatial 
positions of overdense regions (selected manually) to their density weighted positions in velocity 
space. 



3.1.1. Collapsing regions of moderate overdensity 

Figure 4a shows a LOS producing an intermediate column density Hi profile. The line is 

saturated but not yet damped. In LOS passing very close to at least one PGC generally a single 
sharp density peak dominates the absorption line formation (see the following pictures), but more 
often the LOS passes through regions of smaller overdensities just about to collapse into a single 
object. The peculiar velocity field displays a typical infall pattern (redshift with respect to the 
Hubble flow for gas falling in from the front, blueshift for gas falling in from behind, a jump at the 
location of the density maximum). 

Civ and Ovi are the only metals which can be detected. The absorption features are similar 
in appearance to the "partial Lya limit systems" studied by Cowie et al. 1995. In this particular 
case the individual density peaks have converged in velocity space to form a caustic, but the 
enhanced density in physical space is still much more important for producing the absorption line 
(see also Miralda-Escude et al. 1996). In principle velocity caustics may reduce the contribution of 
the Hubble expansion to the total velocity width of absorption lines. We find, however, that this 
rarely occurs. In most cases infall velocities overcompensate the Hubble flow. The ordering of two 
absorption lines in velocity space can even be reversed relative to their spatial positions. 
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3.1.2. Large velocity spread and chance alignments with filaments 



Absorption caused by a chance alignment between the LOS and a filamentary structure is 
seen in figure 4b. The filaments generally lie at the boundaries of underdense regions ("voids") 
and expand with velocities in excess the Hubble flow. This is apparent from the divergent dotted 
lines and the ramp-like increase of the peculiar velocity as a function of real-space distance. Such 
structures rather than deep potential wells are likely to be responsible for the largest observed 
velocity splittings of Civ systems whereas moderate velocity splittings (up to 200 to 300kms~^) 
are often due to a complex temperature distribution or ionization structure internal to individual 
PGCs (see also paper I). 



3.1.3. '^Damped Lya systems 



Figure ^ shows two "damped Lyman a systems" . Currently we are limited to qualitative 
statements about the spatial extent and ionization state of damped Lyman a systems as neither 
self-shielding nor energy /momentum input by star-formation or AGN have been taken into 
account. Nevertheless, in agreement with Katz et al. (1996) we find that for small impact 
parameters column densities becomes large enough to produce damped Lya absorption lines. 
Damped systems are usually dominated by one large density maximum. The density peak 
corresponds to a local temperature minimum due to cooling of the dense gas at the center and 
is surrounded by a shell of hotter infalling shock-heated gas. This leads to the characteristic 
double- hump structure in the temperature diagram of figure 5b. The absorption lines for the 
various ions of the system shown in figure 5a appear rather similar to each other in terms of 
position and line shape, with the exception of Ovi which arises at much larger radii. In general, 
the maxima of the high ionization species for these highly optically thick systems need not coincide 
with the centers of the Hi lines and those of lower ionization species. Figure 5b illustrates a 
situation where the component structure has a very complex origin. While the strongest Cii and 
Siiv components coincide with the center of the damped Lya line, the Civ and Nv positions are 
far off. This indicates differential motion between the high and low ionization regions. 
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Fig. 4. — left: spectrum of a LOS through a coUapsing region of moderate overdensity producing 
intermediate column density Hi and weak Civ and Ovi absorption, right: LOS along a filament 
expanding faster than the Hubble flow. The panels give (from top to bottom) the spectra of the six 
ions offset vertically by 0.5, the total baryon density, the temperature and the peculiar velocity (see 
section 3.1). In the case of Civ, Siiv, Nv, and Ovi only the stronger component of each doublet is 
shown. 
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Fig. 5. — "damped Lya systems" . 
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Fig. 6. — Left: mean projected column density as a function of overdensity (pbar)LOs/Pbar (column 
density-weighted along the line-of-sight) for the simulation box shown in figure 1. Right: column 
density distribution (f(N) N) at z = 3 for six different species in the simulation. Crosses show the 
observed column density distribution (Petitjean et al. 1993). The column density normalization is 
described in section 4.1. 



3.2. The absorption properties of individual protogalactic clumps 



The solid curve in figure 6a shows the good correlation between the mean log A^(Hi) and the 
baryonic overdensity in a typical simulation box. Also shown is the mean log N for our canonical 
set of ionic species. As expected the strength relative to Hi varies considerably from species to 
species. Cii and Siiv are strong at high Hi column densities/baryonic over-densities, while Civ 
dominates at intermediate column densities and Ovi probes the low density regime. We have 
also plotted a simple self-shielding correction for large Hi column densities. The correction was 
calculated with CLOUDY specifying the total column density and mean density along the LOS. 

Figure |^ shows how the spectral features change with distance from a fully collapsed clump. 
The plot consists of a mosaic of LOS separated by 3 kpc from each other in the x and y directions 
on the sky. The center of the clump is close to the top right corner. The Hi Lya line exhibits 
damping wings within the central ~ 6 kpc. 

To quantify the spatial coherence of the absorption properties in more detail we have 
investigated a set of random LOS close to a typical PGC {M ~ 1.2 x IO^Mq) in the simulation. 
Figure ^ shows the column densities as a function of impact parameter for a set of randomly 
oriented LOS. A sharp drop from logA^(Hi) > 20 to logA^(Hi) ~ 15 occurs within the first five to 



15 



ten kpc. At Hi column densities above log 17 self-shielding should become important. This 
will lead to an even steeper rise toward smaller radii. After the rapid drop the Hi column density 
decreases very gradually, still exceeding logA^(Hi) 14 at 100 kpc. The typical radius of the 
damped region of a protogalactic clump in our simulation taking self-shielding into account should 
be about five kpc. However, as pointed out previously this value might increase once feedback 
processes are taken into account. 

Currently a 10m telescope can detect Hi column densities of order 10^^ cm~^ which means 
that the detectable Hi radius extends far beyond 100 kpc. The largest detectable metal absorption 
cross-section is subtended by Civ and probably Ovi (second panel on the left); at a detection 
threshold of A^(Civ)~ 10^^ cm~^ (realistic for very high signal-to-noise data), the radius is ~ 30 
kpc; Siiv at the same significance level (taking the larger oscillator strength into account) would 
be detectable out to about 12 kpc. The decrease of the higher ions Civ and Nv at the lowest 
radii is due to the increase in density; it should be even more pronounced if radiation transfer 
were implemented. Figure |9| gives the ratios of the LOS integrated column densities for several 
ions. In particular, the third panel on the left shows once more the importance of the Ovi ion, 
which should be detectable at the largest radii with similar or even higher significance than the 
Civ lines. Thus, in spite of the difficulty of observing Ovi in the Lya forest the Ovi/Civ column 
density ratio should be the best measure of the impact parameter. 



4. Observational tests 

4.1. Column density distributions as a test of the gas distribution and metallicities 

Before discussing the overall column density distribution of Hi and the metal ionic species 
in our simulations we stress again that our simulation boxes were not chosen randomly (see 
section 2.2). At our redshift of interest {z = 3) the mean baryonic density of the simulation boxes 
containing low (high) circular velocities halos is a factor 1.27 (2.97) larger than the assumed mean 
baryonic density of the universe for the inner 700 (proper) kpc cube. Furthermore, Hi column 
densities projected across the simulation box are generally larger than a few times lO^^cm"^. 
Below we will concentrate on properties of the four boxes containing low circular velocities halos 
at z = 3 for which the overdensity is moderate and which should be closest to a fair sample. 
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Fig. 7. — Mosaic of 6x6 lines-of-sight near a collapsed clump, separated by 3 kpc in each direction 
on the sky. The center of the collapsed object is close to the top right corner. The panels give the 
spectra of the six ions offset vertically by 0.5 for clarity. In the case of Civ, Siiv, Nv, and Ovi only 
the stronger line of each doublet is shown. 
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Fig. 8. — Integrated column densities along random lines-of-sight for the six ions as a function 
impact parameter from the center of a protogalactic clump with 1.2 x 10^ M0. 
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Fig. 9. — Ratios of the integrated column density along random lines-of-sight for several ions as a 
function of impact parameter from the center of a protogalactic clump with 1.2 x 10^ M©. 
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Fig. 10. — Ratios of the integrated column densities along random lines-of-sight for several ionic 
species as a function of Hi or Civ column density. 
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Column density distributions/ratios of ionic species with different ionization potentials are 
an important diagnostic of the ionization parameter, the spectral shape of the ionization field, 
the relative element abundances and the temperature of the gas (e.g. Bergeron &: Stasinska 
1986, Chaffee et al. 1986, Bechtold, Green & York 1987, Donahue & Shuh 1991, Steidel 1990, 
Viegas &: Gruenwald 1991, Ferland 1993, Savaglio et al. 1996, Songaila & Cowie 1996). In our 
case the distributions of densities and temperatures are predicted by the simulation (only weakly 
dependent on the input radiation field) so the remaining adjustable parameters are the elemental 
abundances and the properties of the ionizing radiation. 



4.1.1. Column density distribution functions 



We define the differential column density distribution as usual 

d^AA 

where M is the number of lines and dX = (1 + z^/"^ dz (for (70=0. 5). The curves in figure 6b 
give f{N)N at z = 3.07 for the four simulation boxes shifted in column density by 0.3 dex, 
whereas the crosses show the observed distribution (Petitjean et al. 1993). The data on the metal 
ions are again for a constant metallicity of [Z/H]=— 2. The self-shielding correction for large Hi 
column densities mentioned in section 3.2 is also shown. The shape of the observed distribution 
is matched reasonable well by the simulations once the self-shielding correction is applied. As 
pointed out e.g. by Hernquist et al. (1996) and Miralda-Escude et al. (1996) the Hi column density 
should scale with the baryon fraction, the strength of the ionizing background and the Hubble 
constant as VL^h? / J . With the parameters of our simulation the applied shift of 0.3 dex in column 
density corresponds to (1^^ /iIq/O.OS)^ /{J-22 ^50) ^ 1-5, similar to the value found by Hernquist et 
al. (1996). Taking into account that the mean baryonic density of the volume used to calculate the 
column density distribution is a factor 1.27 larger then the assumed overall mean baryonic density 
should somewhat lower this value. We do also find a deficit of systems with Hi column densities 
around 10^''cm~'^, but it is smaller than the discrepancy by a factor of ten reported by Katz et 
al. (1996). 



4- 1.2. Column density ratios for random lines of sight 



Figure IC shows column density ratios for various ions, integrated along randomly offset and 
oriented LOS through the box described at the end of section 2.3. The most readily observable 
ratio is Civ/Hi, which in the simulations has a maximum of about —2.2 at A^(Hi)=10^^ cm~^ in 
the log. 
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Fig. 11. — Column density ratios obtained from profile fitting the simulated data. Absolute 
metallicity [Z/H]= —2, solar relative abundances and a power law UV spectrum with a=1.5 and 
J_22 = 3 were assumed. Filled circles denote the ratios of the sums of the column densities of all 
individual components along a LOS detectable at the 4.75 a level. The thin diagonal lines give 
upper limits for the simulated lines in cases where only one ion (Hi in the first panel and Civ in the 
other cases) has been detected. Open squares mark observed values (sec text), open circles cither 
observational upper limits or (bottom left panel) Ovi values undetermined because of blending. 
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This is consistent with the range of values measured by Cowie et al. (1995), but a better match 
can be found with a carbon abundance of [C/H]= —2.5 instead of the originally adopted [C/H]= 
—2. At higher column densities Civ recombines to form Cm and Cii with increasing density and 
shielding. At lower column densities Civ/Hi also declines as carbon becomes more highly ionized 
into Cv. This effect will cause difficulties for searches for carbon enrichment in the Lya forest at 
column densities much lower than the presently accessible limits in Hi (10^^ cm~^), as not only 
the Hi column but also the Civ/Hi ratio declines. 

Results obtained from Voigt profile fitting a subset of the spectra (using the software package 



VPFIT, Carswell et al. 1987) are shown in figure 11. Solid dots are measurements from the 
simulation. For comparison open boxes give observational data for Cii, Civ, Siiv and Nv by 
Songaila & Cowie (1996) and for Ovi from Ranch et al. (in prep.). The circles show observational 
upper limits or indicate undetermined Ovi in cases of severe blending with Lya forest lines. 
The thin diagonal lines indicate the approximate detection limits for absorption features in the 
simulated spectra. Most of the general trends in the observed column density ratios are well 
reproduced. The scatter in the column density ratios between different metal ionic species is 
similar to that found observationally. This suggests that the simulations produce a realistic range 
of ionization conditions. 

There are, however, some interesting discrepancies between simulated and actual data. First 
we note that the scatter in Civ/Hi (panel on top left) for the observations is much larger than 
that in the simulations (where constant metallicity was assumed). We take this to imply that the 
carbon fraction and probably also the absolute metallicity in the outskirts of protogalactic objects 
at z ~ 3 fluctuates over one to two orders of magnitude throughout the column density range 
A^(Hi) 10^^ to 10"*^^ cm~^. It is interesting to note that a similar scatter is seen in [C/Fe] in metal 
poor stars (McWilliam et al. 1995, Timmes, Woosley & Weaver 1995). 

There are also some obvious differences between the absolute levels of the observed and 
simulated mean column density ratios. While the abundance ratios for the high ionization species 
(Civ, Nv, Ovi) in the simulation agree reasonably well with the observations, Siiv and Cii are 
significantly off, with the simulated Siiv/Civ and Cii/Civ ratios lower by a factor ten than the 
observed values. 



4-. 1.3. Matching the observed column density ratios 



Some of these discrepancies can be reduced by an appropriate modification of the relative 
abundances. So far we have assumed solar relative abundances. The generally low abundances in 
present-day galactic halo stars (which may contain a fossil record of the high z gas abundances) 
as well as optically thin Civ systems (Cowie et al. 1995) and high-redshift damped Lya systems 
(Pettini et al. 1994, Lu et al. 1996) indicate that relative metal abundances like those in very 
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metal poor stars (McWilliam et al. 1995) may be more appropriate. Below we will discuss the 
effect of changing the abundance pattern from solar values to: [C] = [N]=0; [Si] = [0] =0.4. The 
brackets denote the difference to solar abundances in dex. 

Some improvement may also be expected from a different intensity or spectral shape of the 



ionizing radiation background. Figure 12 shows the influence of such changes on various column 
density ratios. The solid curve shows the values for the a = 1.5 power law spectrum. The other 
curves are for several different spectral shapes and normalizations of the UV background: a 
power law with intensity increased and decreased by factors of 3; a spectrum taking intergalactic 
absorption and emission into account (Haardt & Madau 1996); and a power law with a step-shaped 
break at 4 Rydberg (flux reduced by a factor 100 at 4 Rydberg, constant flux for higher energies up 
to the point where the flux equals that of the underlying power law) . The effect of the absorption 
edges in the Haardt & Madau spectrum on the ion ratios considered here is generally small for 
the relevant column densities. A power law with slope 1.5 is a good approximation. Changing the 
normalization has noticeable but moderate effects on the Siiv/Civ and Cii/Civ ratios. Otherwise 
the changes are again small. The case of a powerlaw with a 4 Ryd cutoff, however, leads to strong 
departures in the ion/Hi, and in the ion/ion ratios for a given Hi column density. A spectrum of 
the last kind was suggested by Songaila & Cowie (1996) and Savaglio et al. (1996) to account for 
large observed Siiv/Civ ratios at z > 3. These authors consider an increasing Heii opacity or an 
increasing stellar contribution to the UV background towards higher redshifts as possible causes 
for a softer spectrum. 

However, for a fixed Civ column density the power law with 4 Rydberg cutoff actually even 
lowers the Siiv/Civ and Cii/Civ ratios relative to those obtained with the originally adopted a 
= 1.5 power law. Moreover, such a spectrum would imply a factor three to ten smaller carbon 
abundance. Apparently, changing the UV radiation field does not dramatically improve the overall 
agreement with the observed ion ratios. 

The observed and simulated column density ratios are compared in figure |l^ after adopting 
some of the changes discussed above (metallicity reduced by 0.5 dex to [Z/H] =—2.5; metal poor 
abundance pattern; spectrum as suggested by Haardt & Madau). The discrepancy for Siiv and 
Cii has now been significantly reduced: the contribution to the shift in the Siiv/Civ vs. Civ came 
in equal parts from the adjustment in relative and absolute metallicity, whereas the Cii/Civ vs. 
Civ ratio improved mostly because of the decrease in absolute metallicity. Thus the agreement is 
now quite good for Siiv, Cii and Civ. Ovi and Nv are slightly further off but ~ given the large 
observational uncertainties in these two ions - probably consistent with the data. 
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Fig. 12. — Column density ratios for various UV background radiation fields. The solid curve 
are for the a. = 1.5 power law spectrum. The other curves are for: a power law with intensity 
increased and decreased by factors of 3; a spectrum taking intergalactic absorption and emission 
into account (Haardt & Madau 1996); and a power law with a step function like break at 4 Rydberg 
(flux reduced by a factor 100 at 4 Rydberg, constant flux for higher energies up to the point where 
the flux equals that of the underlying power law) . 
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Fig. 13. — Column density ratios obtained from profile fitting the simulated data (filled cireles). 
Compared to figure 11 the absolute metallicity has been changed from [Z/H] = —2 to —2.5, solar 
relative abundances have been replaced by those of metal-poor stars ([C] = [N]=0, [Si] =[0] =0.4) and 
the spectrum suggested by Haardt & Madau (1996) is used instead of a power law with a = 1.5. 
The diagonal lines give Tipper limits for the simulated absorption features in cases where only one 
ion (Hi in the first panel and Civ in the others) has been detected. Open squares mark observed 
values (see text), open circles either observational upper limits or undetermined Ovi . 
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Fig. 14. — Left panel: Total Civ Doppler parameter measured by profile fitting versus density 
weighted Civ thermal Doppler parameter computed from the gas temperature in the simulation. 
Right: Total measured CiV Doppler parameter vs. the Doppler parameter obtained by adding RMS 
bulk velocity dispersion and thermal velocity in quadrature. 

4.2. The Doppler parameter as indicator of gas temperatures and small-scale bulk 

motions 

Temperature and bulk velocity measurements from individual absorption lines are useful 
discriminants for the environment of heavy element absorption systems. The actually observable 
quantity is the Doppler parameter b (=\/2cj) of the absorbing line. The line width results 
from a convolution of the thermal motion and the small scale bulk motion. Measurements of 
the Doppler parameter for ionic species with different atomic weights permit a decomposition 
into the contributions from thermal and non-thermal (bulk) motion. The precise nature of the 
decomposition depends, however, on the velocity distribution of the bulk motion. 

For the simulations the relative importance of the contributions to the line formation can be 
easily studied as we have full knowledge of the density and peculiar velocity field along the LOS. 

Here we proceed as follows: first we compute the column density weighted temperature and RMS 
velocity dispersion in an overdense region selected manually; then we fit the absorption line closest 
to the position of the velocity centroid of the region with a Voigt profile and compare results. 
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Fig. 15. — Doppler parameter - column density (6-logA^) diagrams for CiV (top), Ovi (middle) and 
Hi (bottom). The dotted lines give the spectral velocity resolution. Note the different scales of the 
axes. 
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Fig. 16. — Comparison between the Civ Doppler parameters from the simulation (solid line) and 
the observed distribution from Ranch et al. 1996 (dotted line). The observed distribution has been 
scaled such as to match the total number of lines of the simulated distribution. 



The top panel in figure 14 shows a plot of the "observed" Civ Doppler parameter from the fit 
versus the purely thermal Doppler parameter computed straight from the temperature array of 
the simulation. While the density weighted thermal Civ Doppler parameter of many LOS hovers 
around 7 to 8 kms~^, the "observed" values occur over a wide range. The total Civ Doppler 
parameter is obviously not an unbiased measure of the temperature alone. There are substantial 
varying contributions from bulk motion. It is not obvious how temperature and bulk motion 
should be deconvolved. Making the simplest possible ansatz of adding them in quadrature (which 
would be strictly true if the bulk motion followed a Gaussian velocity distribution) we arrive at 
the results shown in the lower panel: the observed line width does indeed measure the quadratic 
sum of thermal motion and RMS velocity dispersion to a reasonable accuracy. Outliers in the plot 
are due small spurious line components sometimes introduced by the automatic fitting program 
VPFIT to improve the quality of the fit. 

How realistic are the gas motions in the simulation? To compare the simulation to real data 
we have again measured a number of randomly selected simulated absorption lines fitting Voigt 
profiles. The resulting Doppler parameter column density (6-logA^) diagrams for Civ, Ovi, and 
Hi are shown in figure 15. The scatter plots agree very well with those from observational data 
(e.g. RSWB for Civ, Hu et al. 1995 for Hi; there is no comparable information yet for Ovi). The 
presence of a minimum b parameter in both cases and a slight increase of b with have both 
been observed. There is also excellent quantitative agreement with the observed mean b values: 
the mean Hi Doppler parameter in the simulations is (6)sim = 27.3 kms~^ ((^)obs = 28 kms~^ is 
observed). The mean Civ Doppler parameters are {b)sim = 8.6 kms^^ and (6)obs = 9.3 kms~^, 
respectively. 

Figure |l^ shows that the shape of the simulated Civ Doppler parameter distribution matches 



29 



20 40 60 

impact parameter [kpc 



20 40 60 

impact parameter [l<pc] 



Fig. 17. — Strongest Civ Doppler parameters of each absorption complex as a function of impact 
parameter to the nearest PGC (left). Mass-weighted temperatm'e (right). 

the observed distribution very well. A decomposition of the measured Doppler parameters 
into thermal and turbulent contributions gives a mean simulated thermal Doppler parameter 



(6)therm = 7.7 km s ^(see figure 14) with the observed mean thermal Civ Doppler parameter, 
(&)thcrm = 7.2 kms"^). 

Can we observe the infalling motion during gravitational collapse directly? Naively we may 
expect to see substantial line broadening due to infalling gas. However, the Hi density weighted 
velocity dispersion along the LOS (the physical quantity relevant for the line broadening) has 
a median as low as 11.6 km s~^, even though the average radius of the regions used for the 
computation of this quantity was as large as 40 kpc. Looking at the peculiar velocity diagrams in 
figures ^ and ^ with their large velocity gradients of close to 200 km s~^ across the size of a PGC 
the quiescent structure may be somewhat surprising. However, the enhanced density responsible 
for the line formation is strongly peaked on a spatial scale much smaller than the infalling region. 
The line profile samples mostly the high density gas which has gone through the shock front, has 
cooled to temperatures of a few 10^ K, and is more or less at rest. In principle the signature of 
infalling gas is still visible in the form of broad profile wings. In Civ, however, these are almost 
always too weak to be seen as they are due to low column density gas beyond the radius of the 
shock. The situation may be different for Hi absorption lines (cf. Ranch 1996). In Ovi with its 
much larger cross-section and weaker density dependence we probably also have a better chance 
to detect infall. The Ovi Doppler parameter ((6)sim=ll-5 kms~^) is indeed higher than that of 
Civ, and not lower, as one could naively expect for thermally dominated gas. Obviously the Ovi 
absorbing gas is subject to stronger non-thermal motions. 

How do the Doppler parameters vary as a function of impact parameter? The LOS-integrated 
temperature distribution of either Civ or Ovi does not show a radial variation that could be 
used as a measure of the impact parameter. This is due to the fact that some of the absorption 
components are related to different PGCs. Plotting, however, only the Doppler parameters of the 
strongest absorption component in each Civ complex (top panel, figure 0) we obtain a noticeable 
anti-correlation of 6(Civ) with impact parameter from the center of the closest PGC. This seems 
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Fig. 18. — Left: two-point correlation function for the Civ systems detected in a simulation box 
at 2=3.07 (to contain three galaxies with Vc ~ 100kms~^at z = 0). Normalization is by the 
combined TPCF from the observations of three QSOs ((z)=2.78; RWSB). Right: same as above 
for a simulation box containing a single galaxy with Vc ~ 200 kms~^ by z = 0. Note the different 
scales of the axes. 

to reflect the drop in temperature (bottom panel) and the coming to rest of the gas at small 
radii (< 12 kpc in this case). The trend is not obvious for Ovi (not shown here) which arises 
from a much larger region of lower density. It also does not show up in plots of 6(Civ) when all 
components, rather than the strongest in each complex, are considered. 



4.3. The two-point correlation function and the dynamical state of the absorbers on 

large scales 

As discussed earlier the observed velocity spreads among absorption components are difficult 
to reconcile with a dynamical velocity dispersion. In paper I it was argued that occasional large 
velocity splitting of sometimes more than 1000 km s~^ can be explained by chance alignments of 
the LOS with filaments containing several PGCs. The observed two-point correlation function 
(TPCF) can indeed be interpreted in this sense, if there is a strong small scale clustering of Civ 
components on scales of 20kms~^ in addition to the large scale expansion (RSWB). This is in 
good agreement with the earlier result that small scale structure may explain the "supra-thermal" 
widths commonly found for high column density Lya lines (Cowie et al. 1995), as well as the 
increasing clustering of Lya forest lines as one goes to higher column density (Chernomordik 1995, 
Cristiani et al. 1995, Fernandez-Soto et al. 1996). 

The TPCF for the simulation box described earlier (containing three Vc ~ 100 kms~^ galaxies) 
is shown in the upper panel of figure ^ The velocities of the Civ systems were obtained by 
Voigt profile fits to significant continuum depressions. The absolute normalization comes from the 
observed TPCF for three QSOs ((2;)=2.78; RSWB). Similar to the observed TPCF, the correlation 
function of Civ lines in the simulated spectra exhibits a narrow peak at the origin (the lowest 
velocity bin is incomplete) and a long tail out to velocities of 500-600 km s~^ and beyond. The 
overall width is somewhat smaller than that of the observed distribution. Repeating the analysis 
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for another simulation box containing one galaxy with v^. ~ 200 km s we obtain the TPCF shown 



in the lower panel of figure 18 which has much more power on large scales. Considering the rather 
small size of our simulation boxes, such variations from box to box are not surprising. In fact the 
observed TPCF is also subject to a large scatter from QSO to QSO, as some LOS often contain 
only one or a few systems providing large velocity splittings. It seems likely that the observed 
TPCF can be explained by averaging over TPCFs from individual galaxy forming regions. 

As a formal measure of the velocity dispersion of an individual protogalactic region we can 
compute the mean (over all LOS) of the RMS of the density weighted total velocity 



I J{v — vYdn \ 



where v = H{z)r + fpec, and 
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For the standard box, with 3 future galaxies, av = 100 km s ^. Four other boxes developing into 
single ~ 200kms~^ galaxies give cr^ = 142, 134, 111 and 138kms~^. 



4.4. Metal absorbers in emission 



The observational identification of Civ absorbers at redshifts z > 2 with any known type of 
galaxies has proven difficult. Speculating that the absorbers may be related to luminous galaxies 
with an old population of stars Aragon-Salamanca et al. (1994) have obtained K images of QSO 
sight lines with known strong Civ absorbers to a K limiting magnitude of 20.3. A slight excess of 
objects probably at the redshift of the QSO was found but convincing identifications of individual 
absorbers were not possible. 

Here we will briefly consider stellar population synthesis models (Bruzual & Chariot 1993, 
Chariot 1996) to investigate the prospects of detecting the stellar continuum from the PGCs 
responsible for the metal absorption systems in our model (see also Katz 1992, Steinmetz 1996a). 
We assume that the available gas turns into stars on time scales between 10^ and 10^ yr, a time 



span capturing the range from a short burst to extended star formation. In figure |19| the apparent 
brightness (not corrected for Lya and dust absorption) at z = 3.1 is shown for a total gas mass 
of 10^ Mq (typical for a PGC) and three different star formation timescales as a function of the 
redshift where star formation began. 



Figure 19 shows that our model would explain the non-detection of Civ absorbers in K 



by Aragon-Salamanca et al. (1994). Even a very short burst of 10^ Mq reaches only a peak K 
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magnitude of 22. Extended star formation reduces the maximum K magnitude to 25. A gas/star 
mass ^ 10^^ Mq would be necessary to be permanently detectable in K from the ground. 

However, if galaxies start building up at redshifts 2; ~ 3 by merging of smaller objects recent 
star formation should have occurred at these redshifts. The optical passbands are then much more 
promising for a ground-based detection of the stellar continuum. Figure |l^ suggests that with 
HST the stellar continuum should be detectable in all wavebands redward of the Lyman break. 

Spectroscopic identification will only be possible for the most massive objects in the case of 
extended star formation, and for a fraction of bursting PGCs of average mass if the star formation 
timescale is short. Recently Steidel et al. (1996) and Giavalisco et al. (1996) have reported the 
detection of a population of star forming galaxies at comparable redshifts. It is intriguing that 
these objects could be identical to either a bursting fraction of average mass (10^ M0) objects or 
to the high mass end of the population of PGCs, showing extended star formation (Steinmetz 
1996a). 

Finally one should note that confusion with foreground galaxies is a concern when studying 
regions of ongoing galaxy formation as described in this paper. The inner 700 kpc of our simulation 
box (at z=3) shown in figures 1 and 3 correspond to an angular size slightly larger than that 
of the WFPC2 field. Thus images as large as the Hubble deep field may be required to get a 
"complete" picture of the progenitors of an individual large z = galaxy. Superposition effects are 
then obviously quite severe. 



33 




Fig. 19. — Bruzual 8z Chariot models for the apparent brightness (not corrected for Lya and dust 
absorption, Salpeter IMF, [Z/H] = —1.7) at z = 3.1 for a total gas mass of a protogalactic clump 
of 10^ Mq, as a function of the redshift where star formation began. Different Hne styles indicate 
star formation timescales between 10^ and 10^ yr. Different panels are for different Johnson filters. 
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5. Discussion and conclusions 



Gravitationally driven density fluctuations in a universe with hierarchical structure formation 
can explain QSO absorption phenomena at 2; ~ 3 over a wide range of column densities. While 
neutral hydrogen shows a rather tight correlation between column density and total density over a 
density range from IQ-^ to 10"^ cm , other ionic species probe different density and temperature 
regimes in a way specific to each species. The lowest Hi column densities (10^^ to 10"*^^ cm^^) arise 
from large-scale sheet-like density enhancements in the IGM with an overdensity of only a few 
compared to the mean density of the universe. In this diffuse gas (densities around 10~^ cm~^) 
high ionization species are prevalent and Ovi A1031 is often the strongest metal absorption line. 
Towards higher Hi column densities we start probing filaments embedded in the large-scale sheets. 
In these regions low column density ClvAA1548, 1550 lines from infalling gas with densities around 
10~^ cm~^ (overdensities of about 10 to 100) dominate the metal absorption features. Civ remains 
the most easily visible metal ion in the as yet unvirialized regions around the protogalactic clumps 
which are later to merge into present-day galaxies. Still larger Hi column densities occur for 
lines-of-sight approaching the central regions of PGCs. These give rise to Lyman limit systems 
and eventually to damped Lya absorbers. Total densities here exceed 10~^cm~^, and species like 
Cll and SilV become increasingly prominent. At densities above 10~^ cm~^ we have reached the 
virializcd region which is generally optically thick for radiation shortward of one Rydberg. 

Although current simulations cannot precisely constrain the size of the damped region in a 
PGC, the fact that at z ~ 3 more than ten such objects exist in the comoving volume containing 
one galaxy at z = considerably reduces the cross section per object required to explain the 
observed rate of incidence of damped Lya absorbers. If this picture is correct there should be no 
1:1 correspondence between present-day galaxies and a high-redshift damped (or Lyman limit) 
absorber. There is then no need for hypothetical large disks/halos as high redshift progenitors of 
present-day galaxies. 

The analysis of artificial spectra generated from our numerical simulation was carried out 
using the same methods as for observational data. The results are in remarkable quantitative 
agreement with a number of observed properties: 

The predicted shape of the Hi column density distribution shows good agreement with that of 
the observed distribution and a good fit is obtained for (fif, /i5q/0.05)^ /{J-22 ^50) ^ 1-5. Detailed 
information on the strength of accompanying metal absorption has recently become available for 
Hi column densities as small as a few times 10^^ cm""^. We obtain good overall agreement between 
the results from our artiflcial spectra and the observed properties, either with a simple power law 
UV spectrum (a=1.5, J-22 = 3) or with the spectrum proposed by Haardt &; Madau (1996). This 
seems consistent with the lower end of the range of J values measured from the proximity effect 
(e.g. Giallongo ct al. 1996, and rcfs. therein) and suggests a baryon fraction slightly exceeding the 
nucleosynthesis constraint (Walker et al. 1991). 
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The scatter of the column density ratios for Siiv/Civ, Cii/Civ, and Ovi/Civ versus Civ is 
consistent with the observational results, so the range of ionization conditions appears to be well 
captured by the simulations. 

As already appearent from simple photoionization models the observed metal line strength 
corresponds to a mean metallicity [C/H] = —2.5 for the column density range cm . 

A homogeneous metal distribution reproduces the observed ion ratios quite well. This implies 
either that much of the as yet unvirialized gas had been subject to a widespread phase of stellar 
nucleosynthesis well before redshift three, or that metal transport outward from fully collapsed 
regions has been efficient. 

Nevertheless, the observed scatter in [C/H] is larger by a factor three to ten than predicted by 
our numerical simulations where the metals were distributed homogeneously. This may indicate 
that some of the metal enrichment took place in situ with incomplete mixing prior to observation. 
Alternatively there may be a wider spread in physical conditions (e.g. spatial variations of strength 
and spectrum of the UV field) than assumed by the simulations. 

The observed column density ratios are matched significantly better if we use relative 
abundances appropriate for metal-poor stars (we used [C]=[N]=0 and [Si] = [0] =0.4). Damped 
Lya systems at redshifts > 3 show similar low metal abundances (—2.5 < [Z/H] < —2.0) and 
relative metallicities (Lu et al. 1996, and refs. therein). This is consistent with the idea that 
metal absorption systems at high redshift contain a record of early nucleosynthesis dominated by 
supernovae of type II. 

The observed distribution of Doppler parameters and the relative contributions to the line 
width from thermal and non-thermal motion are well reproduced by the simulations. Obviously, 
shock heating is a second important heating agent (in addition to photoionization heating) for 
regions of the universe with overdensities between ten and a few hundred. Inspite of the large 
peculiar velocities (~ lOOkms"^) of the infalling gas Civ absorption lines are typically only ~ 
8 to 10 kms~^ wide. This is, because the Civ optical depth arises mostly in narrow post-shock 
regions where the shocked gas has already come to rest and is cooling rapidly. The contribution of 
bulk motions to the Doppler parameters of Ovi and Hi are larger because much of the absorption 
arises at larger impact parameters, where infall of gas is a more important broadening agent. 
In this particular model of structure formation, the CiV Doppler parameter (at z = 3) is, to a 
good approximation, the measure of the quadratic sum of the thermal and the RMS bulk velocity 
dispersion. 

The large scale structure in velocity space, as measured by the Civ TPCF, is consistent 
with the observed Civ TPCF. This supports the hypothesis proposed in paper I that LOS which 
intersect expanding large scale filaments with embedded PGCs contribute significantly to the high 
velocity-tail of the Civ TPCF. The existence of a hypothetical class of abundant deep potential 
wells at high redshift is not required. 
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We have discussed the prospects of detecting the stehar continuum which is expected from 
protogalactic clumps if at least some metal enrichment has occurred in situ. These objects should 
be visible at all optical wavelengths in deep images with the HST and in all optical passbands 
from the ground, at least longward of the Lyman break. We also suggested to interpret the Lyman 
break objects at redshifts around three reported by Steidel et al. (1996) in terms of the high mass 
end of the PGCs causing metal absorption systems or of a bursting fraction of lower mass PGCs. 
Prospects for identifying metal absorption systems at these redshifts are good. The best strategy 
may be a systematic search for Lyman-break objects within a few arcseconds to the line-of-sight to 
bright quasars, together with follow-up spectroscopy to the faintest possible limits. PGCs which 
are progenitors of a particular z = galaxy can be scattered over several hundred kpc at z = 3, 
an area larger than the field size of the Hubble WFPC2 camera, so problems with incompleteness 
and foreground confusion can arise. 

In future, a large database of metal-line ratios as a function of redshift should improve 
constraints on the normalization and spectrum of the UV background and allow us to distinguish 
between possible metal enrichment histories of the IGM. In particular, in the low density regions 
probed by the lowest column density Lyct absorbers, we expect Ovi to be considerably stronger 
than Civ. Thus it may be possible to push mctallicity determinations with Ovi closer to truly 
primordial gas than is possible with Civ, despite the severe problems with identifying Ovi in 
the Lya forest. Another interesting case for further study is Siiv. Siiv should rapidly decrease 
towards low densities and a detection in low Hi column density systems would indicate that the 
normalization and/or the spectrum of the UV background differ significantly from those we have 
adopted. A detection of the stellar continuum expected to be associated with metal absorption 
systems would allow us to investigate the radial density and metallicity profile of the PGCs. 

In summary, we have found that high-resolution hydrodynamical simulations of galaxy 
forming regions can substantially aid the interpretation of the observed properties of metal 
absorption systems. Currently these simulations are best suited for studying regions with 
overdensities from ten to a few hundred. Such regions are optically thin to ionizing radiation. 
Substantial simplifications are present in the current work. However, the good and sometimes 
excellent agreement of the present model with observation gives reason to believe that we may 
have correctly identified the mechanism underlying many of the metal absorption systems at 
high redshift; conversely, we may take the results presented here as an argument in favor of a 
hierarchical galaxy formation scenario. 
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